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Abstract 



We consider the asymptotic behavior of the second mixed moment of the char- 
acteristic polynomials of the ID Gaussian band matrices, i.e. of the hermitian ma- 
trices H n with independent Gaussian entries such that {HijHik) = SikSjlJij, where 
J = (-W 2 A + I) -1 . Assuming that W 2 = n 1+e , < 6 < 1, we show that this 
asymptotic behavior (as n — > oo) in the bulk of the spectrum coincides with those 
for the Gaussian Unitary Ensemble. 

m 

1 Introduction 

in 
oo 

The hermitian Gaussian ID random band matrices (RBM) are hermitian (2n+l) x (2n+l) 
matrices H n (we numerate indices of entries from — n to n) whose entries Hjj. are random 
^ . Gaussian variables with mean zero such that 

(N . 

(HijHik) = 5 ik 5jiJij, (1.1) 
where Jij is a symmetric function which is small for large \i — j\ and 



E J *i = L 

i=—n 

In this paper we consider the especially convenient choice of J^, which is given by the 
lattice Green's function 

J ij= (-W 2 A + l)-\ (1.2) 
where A is the discrete Laplacian on [— n, n], i.e. 

(-Af)j = I ~ f r + 2 /> ~ 3 * (1-3) 
v 13 I -fj-i + fj ~ fj+i, J=-n,n. y ' 

Note that for this choice of J we have p~ W~ x exp{ — \i —j\/W}, and so the variance of 
the matrix elements is exponentially small when ^> W. Hence W can be considered 

as the width of the band. 
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The probability law of the RBM H n can be written in the form 

dH dH* \ H a\ 2 n dH h 2 . 

PMHn)= n ^« c -^ n ^ e -*, (L4) 

-n<i<j<n A7[J% 3 i= _ n V**J n 

We assume also that W 2 = n 1+e , < 6 < 1, i.e. the width of the band tends to 
infinity faster than y/n, when n — > oo. 

Let . . . , be the eigenvalues of H n . Define their Normalized Counting Measure 
(NCM) as 

N n (A) = eA,j = l,...,n}/n, N n (R) = 1, (1.5) 

where A is an arbitrary interval of the real axis. For the RBM it was shown in [15] that 
N n converges weakly to a non-random measure N which is called the limiting NCM of 
the ensemble. The measure iV is absolutely continuous and its density p is given by the 
well-known Wigner semicircle law: 

p(A) = J_V73p, AG [-2,2]. (1.6) 

Characteristic polynomials of random matrices have been actively studied in the last 
years (see e.g. [1, 4, 5, 7, 9, 10, 12, 14, 16, 17, 18, 19]). The interest to this topic is 
stimulated by its connections to the number theory, quantum chaos, integrable systems, 
combinatorics, representation theory and others. 

The second correlation function (or the second mixed moment) of the characteristic 
polynomials is 

F 2 (A) = J det(A x - # n )det(A 2 - H n )P n (dH n ), (1.7) 

where P n (dH n ) is defined in (1.4), and A = diag{Ai, A 2 } are real or complex parameters 
that may depend on n. 

We are interested in the asymptotic behavior of (1.7) for matrices (1.1) - (1.4) as 
n, W -> oo W 2 = n l+e , < 9 < 1, and for 

A, = A + ^-, j = 1,2, 
np(A ) 

where Ao € (—2, 2), p is defined in (1.6), and £ = diag {£i, £2} are real parameters varying 
in the interval [-M, M] C R. 
Set also 

The main result of the paper is the following theorem: 

Theorem 1. Consider the random matrices (1.1) - (1-4)- Define the second mixed mo- 
ment F 2 of the characteristic polynomials as in (1-7). Then we have 

lim D 2 -WA + f/MAo))) = sin ^-y) , (1.9) 

where p(A) and D 2 are defined in (1.6) and (1.8), A = diag{Ao,Ao}, Ao G (—2,2), 
i = diag {6,6}, £1,6 G[-M,M]. 
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The theorem shows that the above limits for the second mixed moment of the char- 
acteristic polynomials for the ID Gaussian random band matrices (with W 2 = n 1+9 , 
< 9 < 1) coincide with those for the Gaussian Unitary Ensemble, i.e. the local behavior 
of the second mixed moment in the bulk of the spectrum is universal. 

The paper is organized as follows. In Section 2 we obtain a convenient asymptotic 
integral representation for F 2 , using the integration over the Grassmann variables and 
the Harish Chandra/Itzykson-Zuber formula for integrals over the unitary group. The 
method is a generalization of that of [3, 4] and is an analog of the method of [16, 17], 
where the hermitian Wigner and general sample covariance matrices were considered. 
In Section 3 we prove Theorem 1, applying the steepest descent method to the integral 
representation. Appendix is devoted to the proofs of the auxiliary statements. 

We denote by C, C±, K, etc. various n-independent quantities below, which can be 
different in different formulas. Integrals without limits denote the integration (or the 
multiple integration) over the whole real axis. 

2 Integral representation 

In this section we obtain the integral representation for F 2 of (1.7) by using integration 
over the Grassmann variables. This method allows us to obtain the integral representation 
of the product of the characteristic polynomials which is very useful for the averaging 
because it looks like the Gaussian-type integral (see the formula (2.7) below). After 
averaging over the probability measure we can integrate over the Grassmann variables to 
obtain the usual asymptotic integral representation which can be studied by the steepest 
descent method. 

The integration over the Grassmann variables was introduced by Berezin and widely 
used in the physics literature (see e.g. [2] and [6]). For the reader convenience we give a 
brief outline of the techniques. 

2.1 Grassmann integration 

Let us consider two sets of formal variables {V'jljLi, {^jYj^i, which satisfy the anticom- 
mutation conditions 

ipjipk + i^ktpj = i>ji>k + i>ki>j = i>ji> k + ^k^j = o, j,k = i,...,n. (2.1) 

These two sets of variables {^-}™ =1 and {■0 J }" =1 generate the Grassmann algebra 21. 
Taking into account that ip 2 = 0, we have that all elements of 21 are polynomials of 
{■0j}j=i and {i)j}j=i- We can also define functions of the Grassmann variables. Let x be 
an element of 21, i.e. 

n 

x = a + Y^Mi + h $i) + J2( a J,k^k + bj^^k + c jt #ij$ k ) + .... (2.2) 

3=1 j^k 

For any analytical function / we mean by f(x) the element of 21 obtained by substituting 
X — a in the Taylor series of / at the point a. Since x is a polynomial of }™ = i, {^j}" = i 
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of the form (2.2), according to (2.1) there exists such I that (x ~ a Y = 0> an d hence the 
series terminates after a finite number of terms and so f(x) £ 21. 
For example, we have 

exp{6^j} = 1 + b^jipj + = 1 + brj}^, 

exp{a 11 ^ 1 ^l + 012^1^2 + 021^2^1 + «22^2^2} = 1 + OllV'lV'l 

+ 012^1^2 + 021^2^1 + ^22^2^2 + (ailV'lV'l + 012^1^2 (2.3) 

+ 021^2^1 + ^Vv^) 2 / 2 = 1 + On^^i + ayity-^i + a 2 iip 2 ' l l J i 

+ 022^2^2 + (ail022 ~ 012021)^1^1^2^2- 

Following Berezin [2], we define the operation of integration with respect to the anticom- 
muting variables in a formal way: 

J dipj = J dfy = 0, J ijjd^jj = J i^jd^j = 1. (2.4) 

This definition can be extended on the general element of 21 by the linearity. A multiple 
integral is defined to be a repeated integral. The "differentials" d ipj and d ip k anticommute 
with each other and with the variables ipj and ip k . 
Thus, if 

k 

f(m, ■■■,Vk)=Po + ^2 PhVh + Yl P3132V31V32 + ■■■+ pi,2,...,kVi ■■■m, 

jl=l jl<32 

where rji, . . . , rjk are some elements from the sets {■?/>.,• }™ =1 , {^}" = i, then 

y • • • , %)^% ...drn= pi >2 ,...,fc. (2.5) 

Let A be an ordinary hermitian matrix. The following Gaussian integral is well-known 

dtftzid^z-j 1 



/ 



exp 

j,fc=i " j=i 



One of the important formulas of the Grassmann variables theory is the analog of this 
formula for the Grassmann algebra (see [2]): 

/n n 
exp { ~ A i,^k }Y[d^ j d ^ = det A. (2.7) 

A In 1 A 1 



j,k=l j=l 



For n = 1 and n = 2 this formula follows immediately from (2.3) and (2.5). 
Besides, we have 

/q n n 

n^P^p ex P { ~ Yl A j- ! ' '.i' ' 1 '}]!' 1 ^/ 1 l 'i = detA iu-Usu..,s q , (2.8) 

^ 1 ~ 7„ 1 A 1 



P=l j,k=l j=l 



where Ai limm j . slimmiS is a (n — g) x (n — q) minor of the matrix A without lines li, . . . , l q and 
columns Si, . . . ,s q . 
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2.2 Formula for F 2 

Using (2.7) we obtain 



2 -2 n 



F 2 (A) = E e J] J] d^Vv 

^ J r=l g=-n ' 

2 n 2 

= e{ I exp [ £ J2 ■ + (2 ' 9) 

_ _\™ 2 _ "| 2 n _ ^ 

+i^Hj,k ■ {tpji^ki - ^M^ji)) + E Hjj ■ \ EI II difj qr dip qr >, 

' j=—n 1=1 J r=l q=—n ) 



where {ipji}, j = —n, . . . ,n, I = 1,2 are the Grassmann variables (2n + 1 variables for 
each determinant in (1.7)), and E{. . .} is an expectation with respect to the measure 
(1.4). Integrating over the measure (1.4) we get 



2 n 



*2(A) = / I! II dV^exp{ -J2Kj2^ P s} (2.10) 

r=l q=—n s=l p=l 

n j 

X expj Jj,k(^jl1pkl + ^Mi^kl^jl + 'tfkl&j*) + Yl -jf^j^n + ^j2^2) 2 }- 



j<k j=—n 

Use the well-known Hubbard-Stratonovich transform: 



/ ex p { - 1 e - < E(^i» (^) } n tx* 

j,k j j=-n 

= (27r 2 ) 2n+1 det 2 J • exp {i Jj^jiipki + V^Wfe^i + feM}, 



where Xj is hermitian 2x2 matrix and 



dX^- = <mX 12 d3X 12 dX u dX 22 . (2.11) 

Substituting this and expression (1.2) for J^ 1 to (2.10) and using (2.7) to integrate over 
the Grassmann variables we obtain 

F 2 (0 = -( 27 r 2 )- 2 -Met- 2 J / exp { - \ £ Tr (X, - X^) 2 - 1 £ Tr X 2 } 

^ 2 j=-n+l j=-n 

n n 

x Yl det (Xj - iA - ii/np(\ )) JJ dX^- (2.12) 



j=-n j=-n 

^2 



= - (2vr 2 )---Met- 2 J | exp { - \ £ Tr (X, - X^) 2 } 



Hi 

r l A „ 

X 



j=-n+l 



1 ' A '£ 2 

exp { _ i ^ Tr (a-. + f + < ) } n (a-. - iA 0/ 2) n 



np(A ) 

j=—n 3=—n J=—n 



Let us change the variables to Xj = U*AjUj, where Uj is a unitary matrix and 
Aj = diag {cij, bj}, j = — n, . . . , n. Then dXj of (2.11) becomes (see e.g. [13], Section 3.3) 

7T 

-(a, - bj) dajdbjd^Uj), 

where d/i(Uj) is the normalized to unity Haar measure on the unitary group U(2). Thus 
we get 

H A ° + ^j) = - / exp { - ? J?^ {u;m ~ U ^ AM *) 

1 n • A 2 . n 

J=— n j=—n 



n n 

X 

fc=— n 



JJ (a fc -iA /2)(6 fc -iAo/2) II (ai - k) 2 dadbdfi(U q ), 



where 



cia= ]"[ daj, rffe= ] [ c% (2.13) 



j=-ra J=-n 



Note that we can change "angle variables" Uj to Vj = UjU*_ 1: j = —n + 1, . . . , n (i.e. the 
new variables are £/_„, V- n+i , V- n+2 , ■ ■ ■ , V n ). Then we have 



[7(2) 3—n+l 

xexp{-| E Tr(^ + ^) 2 + ^ E T* (£UP,-)% (EUPiK} 

^ j=~n j=—n J 

x II [(a k -b k ) 2 (a k -i\ /2)(b k -i\ /2))dfi(U- n )dadb f] ^(V^), 

k=—n ^ ' p=— n+1 

(2.14) 

where 

p fc = n ^- ( 2 - is ) 



S = — 71+1 



3 Saddle-point analysis 

3.1 Sketch of the proof 

We can rewrite (2.14) as 

-2 T /• , U/2 



xexp{ - E (/(%) + f(bj)) + t^T) E Tr(f/_„P,) Aj{U_ n Pj)i) 

n n 

x [I (a l -b l ) 2 dfj.(U-. n )dadb n 



where 

f(x) = (x + i\ /2) 2 /2 - log(x - iAo/2). (3.1) 

Note that 

d -f(x) = x + t\ /2- ' 



dx x — «Ao/2' 

Hence, the expected saddle-points are 



o± = ±^\3 = ±7r p(A ). (3.2) 

Let E be the integral of (2.14) over the domain Q$ which is the union of 5- neighbor- 
hoods of the points (a + ,a + ), (a + ,a_), (a_,a + ), (a_,a_), where a± = (a±, ...,a±) G 
R 2n+1 , and let E c be the integral over the complement of Cls. 

We are going to prove that 

|E C | < Ce- wl - e I , (3.3) 

where e is sufficiently small and 

J ~ C W" 8n - 2 (27c) 2n+1 2 2n e^ 2n+1 ^det- 1 (-A + ^ , (3.4) 

where 70 is a real positive constant, and Co = 1/2 — Aq/4. 
Since below we will prove that 

E ~ de^^Io, (3.5) 
we obtain then (recall that W 2 = n l+e ) 



V I o-W 1 -"- T 

■ ' < c - - < C e~ wl ~ £/2 



|E + E c | - |e- wl " £ - Cie°W w )|/ 
and thus 

f '( a » + 4j) = -^ s(i+o<i)) - 

As we will prove below, we can take <5 = W^ _K with sufficiently small k. 

The next step is the calculation of E. We take the W^~ K -neighborhood of the one of 
the points (a+,a+), (a + ,a_), (a_,a+), (a_,a_) (for example (a+,a_)). The idea is to 
expand 

/(x) = c ± (x - a±) 2 + s 3 (:r - a±) 3 + . . . = c±(x - a±) 2 + tp±(x - a±), 

C± - 1 "T ± 1"-V 1_ T' 

then to leave only the quadratic form in the exponent in (2.14), and then integrate over 
the unitary groups. After the integration and changing of the variables aj — > a± + dj/W, 
bj — > a± + bj/W , we will face with a problem to study a complex value Gaussian (2n+ 1) - 
dimensional distribution 

Y n n 
/i 7 (x) = exp { - - J2 ( x i ~ x i-i¥ x2 A ( 3 - 6 ) 

j=-n+l j=—n 
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with 7 e C, 9^7 > (in our case 7 = c + or c_). We take 5 = W K and define 



/aw " paw " 

(•••) •/'-(•'•) I f dx q , Z Sn = / //.(./•) TT r/.r, ; (3.7) 

^ ,=-„ ■'-W q =-n 

/n „ n 

(•••)• a*7(z) n ^7 = / n dx <i- 



n 

{ ^^(^/^)}-l) o = o(l). (3.8) 



q=—n q=—n 

Moreover, let (. . .)* (and also (. . .) ,*) be (3.7) with n^{x) instead of // 7 (x). 

To leave only the quadratic form in the exponent in (2.14), we have to prove that 

n 

cxp ' 

j=-n 

This can be done by using three ideas: (1) we can replace the integral over the neighbor- 
hood by the integral over all real axis with an error which we can control; (2) using the 
Wick's theorem we can prove that 

l(4 1 ---4 i >l<(4 1 ---4 i >* (3-9) 

and thus we can estimate the averaging of some function over the complex measure by 
the averaging of the "changed" function over the positive one (we take 5Rc± instead of 
c±); (3) using the Wick's theorem we can prove (3.8) for the positive measure. 
These ideas also help to prove that the main contribution to the integral of 

• 3 3 

exp{— TV(C/_ n J] V s Y(L + A s /W)(U- n J] V s )£} 

^ °> s=-n+l s=-n+l 

gives the term 

exp {^k) Tr(f/ - n vsYHu- n fi v.)i}. 

PV u; s=-n+l s=-n+l 

Here L = diag{a + , a_}, A s = L + A s /W. 

The last step is to prove that in the functions above we can replace all V s by I 2 . 

3.2 Proof of (3.3) 

According to (2.14) we have 

/1 n n 
exp I - -ft ^ Tr (Aj + iA /2) 2 + 3? lo S det (^i - iA o/ 2 )} 

3=~n j=~n 

xex p{-— E Tr W-Vi) 2 } 

j=-n+l 

n n 

x J" [ (oj - k) 2 dfi(U- n ) dadb |] d/j,(V p ). 

l=—n p=— n+1 

The integral over the unitary group U (2) can be computed using the well-known Harish 
Chandra/Itsykson-Zuber formula (see e.g. [13], Appendix 5) 
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Proposition 1. Let A be the normal p x p matrix with distinct eigenvalues {aj}f =1 and 
B = diag{bi, . . . , b p }. Then 

J J expj - ^Tr (A -U*BU) 2 } A 2 (B)f(B)dUdB 

u(p) n 

= (n^' ! ) Jew{-\^(a J -b 3 ) 2 }^f(b 1 ,...,b p )dB, (3.10) 

where f(B) is any symmetric function of {bj}^ =1 in the symmetric domain Q, dB = 
p 

]~[ dbj and A(A), A(B) are the Vandermonde determinants for the eigenvalues {aj}f =1; 
j i 

{hY i= i of A andB. 

We get (recall that Aj = diag {aj, bj}, j = —n, . . . , n) 

Q C j=-n+l j=- n +l 

n 

x exp | - ^2 (f*( a j) + f*( b j))} \ { a ~n - b- n )(a n - b n )\ dadb 

j=-n 

» n n 

=22 ^- 8 „- V ( 2 n +1)co / exp |_l (a.-a^Y- 1 - J2 (b 3 -b 3 -if} 

WU9 j=-n+l j=-n+l 

o 

n 

x exp { - (U a il W ) + Mbj/W))} |(a_ n - b_ n )(a n - b n )\ dadb, (3.11) 

j=—n 

where 

Ux) := ((x + t\ /2) 2 /2 - log (x - i\ /2)) - c , c := ± - ^. (3.12) 
We need 

Lemma 1. The function f*(x) for x G R attains its minimum at x = a±, where a± is 
defined in (3.2). Moreover, f*(a±) = and ifx^ U$(a±) := (a± — 5, a± + 5) for sufficiently 
small 5, then 

U(x) > C5 2 . (3.13) 

In addition, we have for x G (— oo, 5) 

f*{x)>a{x-a_) 2 , (3.14) 

where a is some positive constant. Similar inequality holds for x G (— 5, +oo) (with a + 
instead of a_). 
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Consider a_ n , . . . , a n . First note that the integral over the configurations of {%}"=_„ 
such that I a j ex i j ^ | > Wn £ with sufficiently small e satisfies bound (3.3). Indeed, 



j=-n+l 



- a, ) 2 > CW 2 n 2 = C// 1 " 



Besides, Lemma 1 yields 



f*(x) > a (x — a-) 2 , x < a_, 
f*(x) > a (x - a + ) 2 , x < a + , 



and hence the integral in (3.11) over Y\^ = _ n da q can be bounded by 
C x W 2n + l exp{-Cn 1+e - 2£ } < exp{-Cn 1+e - 2£ /2} 

for any e < 9/2. 

Thus we need to study the configurations such that \dj — a_,-_i| < Wn~~ £ , 
j = —n + 1, ... ,n. Without loss of generality let a_ n < 0. Let li be the first number 
such that a/j > Consider the nearest to l± indexes p± < l\ and q\ > l± such that 

a Pl < 0, a qi < 0. The sequence a Pl+ \, . . . , a gi _i we will call "the peak". Remove from the 

n 

sum ( a j ~ a j~i) 2 the terms (a Pl+1 — a Pl ) 2 and (a qi — a gi _x) 2 (the integral becomes 

j=-n+l 

larger). Then take the first number l 2 > qi such that a/ 2 > SW and the nearest to Z 2 
indexes p 2 < h and q 2 > h such that a P2 < 0, a g2 < and again remove the terms 
(a P2+1 — a P2 ) 2 and (a q2 — a g2 _i) 2 , and so on (the last "peak" can be from a Pj+1 to a_ n ). 
Assume that we obtain k of such "peaks" . 

Consider one of them. Let it has m -\- 1 positive numbers Op r ^_x, . . . , a Pr j rm j r \ — a qr —\. 
Since \a,j — < Wn~ £ , we have m > n £ 5 and taking into account that a Pr < 0, we 
have a Pr+ i/W - a + > 5, \a ir - a Pr+ i| > 5W/2. Set 

m m 

^\x) = exp { - - J>, - x^) 2 - J- J> 2 }. (3.15) 

3=2 3=1 

Lemma 2. For any 7 G C, > we have 

■= j H { ™\x)f[dx q = (27rr/ 2 (^sinh^)~ 1/2 (l + o(l)) (3.16) 



9=1 

Moreover, if we set 



G (m) (7)= ("A + -^) \ (3.17) 



then 



G'r'h) <^tanh-^(l +o( l)) (3.18) 
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(2) 



I y (m) 7 (m)\ 
l Z 7 ~ L & n I 

I ry(m) I 



'H|-l 



max |xi|><5W 



9=1 



-C 2 <5 2 W 



W ->■ oo 



where m > CW , 5 = W K for sufficiently small k. 
In addition, for any m 



| Z (m)|-l 



") n dx « 



9=1 



(3) Let m > C{W , k < Cm/W , S = {i±, . . . , i s } C {1, . . . , m}, and hi = 3fc, where 



h e {3, . . . , k}. Then 



1^7 I 



jJ( Xj /Wp-^\x)l[dx, 



max |xi|>5W 



9=1 



< e 



-d5 2 W 



, m, W — > oo, 



where 5 = W K for sufficiently small k. 

Since a Pr+ i, . . . , a Pr+m+ \ > 0, according to Lemma 1 we can write 

f*(a Pr+s /W) > a (a Pr+s /W - a + ) 2 , s = 1, . . . , m + 1. 

Using this inequality in the r.h.s. of (3.11), we can apply Lemma 2 to the integral over 
da Pr+1 , . . . , da Pr+m+ i. This gets the bound (2n) m / 2 W exp{-m v / 2c^/(2W)-C5 2 W} (recall 
that a Pr+ i/W — a + > 5 and \a [r — a Pr+ i\ > 5W/2). Hence, for integrals over variables 
which make k "peaks" we obtain the bound 

(27r) Emi/2 W k exp{-V2c^J2 m i/( 2W )} exp{-C5 2 Wk}. 

By the same way we can estimate the integral over a qv . . . , a Pl+1 (i.e. over Oj's that lying 
between two "peaks" ) by (2ix) s / 2 W exp{— s\^2a / (2W)} , where s = pi + i — qi + 1. Finally, 
the whole integral with k "peaks" can be bounded by 

(27r) n W 2k+1 exp{-V2an/(2W)} exp{-5 2 Wk}. 

The number of such configurations is smaller than ( 2n 2 ^ 1 ) (since the number of choices of 
the "beginnings" and "ends" of k "peaks" is ( 2n 2 ^ 1 ))- Hence, we get the bound for the 
integral over configurations with k > in (3.11): 



Z- 1 ■ {27i) n exp{-V2c^n/{2W)}J2 



k=l 



2n 



2k 



W 2k+l exp{-5 2 Wk} 



< VWe~ cPw+Cn / w ((l + We~ pw ) n - 1) < e~ c ^ w . 

Therefore, the main contribution to (3.11)is given by the configurations without "peaks". 
For such configurations 

/*(%) > a (% - a+) 2 , j = -n,...,n, 
and thus Lemma 2 proves (3.3). 
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3.3 Calculation of S 
3.3.1 £± and S T 

Consider the ^-neighborhood of the point (a+, a_) with a± of (3.2) and 5 = W^~ K 
We can write for x G Us(a±) 

f(x) = ^(x + i\ /2) 2 - \og{x - i\ Q /2) - r± 



=: c±(x - a±) 2 + (p±(x - a±). 



(3.19) 



where 



c ± = 1 - ^ ± ^° • ^ 1 - ^, r ± = i(a ± + *A /2) 2 - log( a± - *A /2). (3.20) 

It is easy to check that 9ftr± = c , where c is defined in (3.12). 

Doing the change a,j — a + = dj/W, bj — a_ = bj/W in (2.14) we obtain (recall that 
a± = ±7rp(A )) 

E ± =W r-2(2n+l) 2 2n e 2(2n + l)co+< 1 rKi-6) ^ ^ ( a )^_ ( & ) 



M,N<wi-«tf(2) 

E H:(V7(L+^/W)15(L+A J -_i/W)-(L+A J -/W)(L+A J -_i/W)) 

x e 



(3.21) 



- E (¥'+(a fc /VF)+ ¥ p_(6 fc /VF))+ 

x e fc =~ n 



np(A ) 



E (Tr(l/_ n P p )*(L+A p /W) ({/_ n P p )|-TrL|) 



X 



J] (a + - a_ + (5, - 6,)/W) 2 dM^-») II <W)dad&(l + o(l)), 

<7= — 71+1 



l=—n 



where L = diag {a + ,a_}, Aj = diag {dj, bj}, and /i 7 (a) is defined in (3.6). 

We will use bellow the following form of Wick's theorem: for any smooth function / 



( x h f ( x h i ■ ■ ■ i x i P )) ~ ^ ] { x i\ x ij ) (9f (x^ , . . • , Xi p ) / 8. 

3=1 



(3.22) 



The same is valid for (. . .)*, where (. . .), (. . .)* are defined in (3.7). Set 

M = -A+7/W 2 = (2+ 1 /W 2 )I-M, M* = -A+^/W 2 = (2+^/W 2 )I-M, (3.23) 
where 

/ 1 1 ... \ 
1 1 ... 
1 1 ... 



M 



\ 



... 1 1 
... 11 



Then 



(xiXj) = (M (xiXj), = (M„ %. 
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Besides, since all entries of M are positive and 9?7 > 0, 



\(M- 1 )--\ = I V I 

< f < f (^k = (M -i).. 

/c — fc — 

This and (3.22) yield (3.9). 

First step: calculation of ( exp < — Y ( P±{ x j/W) \ ) 

\ j=—n / 

Define 

n 

£„[<?]:= exp { - ^ g(x s /W)} (3.24) 

j=-n 

for any function g : M. — > C. 
We are going to prove that 

|<£„b±]) -l|=o(l), (3.25) 

where y2± are defined in (3.19). 

The key point of the first step is the lemma: 

d 

Lemma 3. Let g = Yl c j x3 ' > c j e ^- Then we have 

J=3 



;^]> 0j# -l|=o(l), n^oo. (3.26) 
Proof. Since e 21 — 1 > x, we have 

n n 

( E nl9])* fi ~ 1 > ( E ^V^o = ( E ^V^)). + 0(1), 
j=-n j=-n 

where in the last equality we use the third assertion of Lemma 2. Using the Wick's 
theorem (3.22) and (M" 1 )^ = CW (see the assertion (1) of Lemma 2), we can write 

{{ Xj /wf l )^o{w- l i 

and hence 

< E 9MW)) m = 0((2n + 1)/PU 2 ) = o(l). 

j=-n 

Let us prove that 

( E »\d)* fi ~ 1 ^ £ m(^]>*,o =► (E n \g]) mfi - 1 < 2e 1>B , (3.27) 
where £i jTl = o(l), as n -> 00. 



13 



Note that if we choose s K > 3 such that 



W~ KSk < W . 



(3.28) 



then for any p > s K /2 and for Xj G (— 5W, SW) 

Y,(^/W) 2p < n/W 2 = o(l), 



and thus if we replace g(x) by g(x) — Cx 2p with any C, then E n [g] will be changed by 
E n [g](l + o(l)) . Since it is easy to see that we can choose C such that cqx 2 /2 + g(x) — Cx 2p 
has only one minimum x = in R, without loss of the generality we can assume that 
c x 2 /2 + g(x) > c x 2 /4:. Moreover, c x 2 /2 + g(x) < c x 2 for x G (—5,5). This and 
assertions (1), (2) of Lemma 2 give 



E n [g}fx C0 (x)dx 



max |xj|>5W 



< 



max |xj|>5W 



fx C0 / 2 (x)dx 



/ E n \g]fjL C0 (x)dx / n 2co {x)dx 

J max \xi\<5W J max\xi\<SW 

because 5 = W~ K with k < 9/2. Thus, 

(En[g])^ = (E n [g])^ + o(l). 

Since for x G R 

< 1 + xe x , 

we can write using the Wick's theorem (3.22) 



< e Cn/W-W&* = o(1)) 



d I 
C l X u 



(En[ 9 ]i -i<E <^/^) • = EE \ T p 

ii ii i=3 



•£n[0] 



<EE-" 1)|Q|(4) '~'~ 2 

ii Z=3 



W 2 



d I I / v Z-l 



i\,%2 1=3 
d 

sEE 

ii 1=4 
d 

+ EE 

ii,«2 i=3 

+ EE 



1^2 



(Z-1)(Z-3)|q|«>2 /x!- 4 



w 4 

(21 - 3)|ci|(a: il a: ia )*(a: 



w l - 



■ E n [g] 



W 4 

|q|(^i^ 2 )* / 



/ Xil n* ( ^11 

\ w l ~ 3 ' y \w 



■E n [g] 



w 4 
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\w l ~ 2 



+ E E 

jl,i2,*3 '=3 



Q I (^1^2) * (^1^3 )* / *^ii J ( X 



w 4 



1 I 



(3.29) 
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Every time having (x™ 1 . . . x™ k E n \g])* we take x i; with the smallest index I and find the 
pair to it according to (3.22). In such a way we make 2ds K steps, where s K is defined in 
(3.28) and d is a degree of the polynomial g. All terms have the form 

G( x hi • ■ • > x i p+ i) 

h,—,i p +i 

where a p+ \, . . . , a p+ i G N are bounded by some absolute constant (since we make the 
finite number of steps), a = a p+ i + . . . + a p+ i, and G(xi 1 , . . . , i p+ i) is the product of the 
expectations of some partition by pairs of x^x k * . . . x*^ with kj > 3, j — 1, . . . ,p and 
kj > 1, j — p + 1, . . . ,p + I (all k/s are bounded by some absolute constant), with some 
bounded positive coefficient and divided by W kl+ '" +kp+l . These partitions we can imagine 
like graphs with vertex i\, . . . , i p+ i, where p + I < 2ds K and the degree of each vertex ij 
is grater than 3 for j < p and is grater than 1 for j — p + 1, . . . ,p + I (these graphs can 
contain loops, i.e. edges like {ij,ij), and multiple edges). Let H be one of such graphs. 
Any (xiXj)* gives (M* 1 )^. Thus, any loop gives factor (M*" 1 )^ = CW(1 + o(l)) (see the 
assertion (1) of Lemma 2). Moreover, according the Cauchy-Schwarz inequality we have 

(M~% < (M-*)f (71C% /2 . 

Thus, we can remove the edge {31,32) from any cycle {31,32, ■ ■ ■ , jr, ji) { r 7^ 1) and change 
it into two "semiloops" (ji,ji), (j2 ? J2) (i-e. we must remember that these "semiloops" 
give {{M-%1 1 / 2 instead of {M'% and one "semiloop" gives the contribution one to the 
degree of the vertex, two "semiloops" in the same vertex give the loop). In this way we 
transform graph H to the collection of the trees H with some loops and "semiloops" 
(the degree of each vertex still the same as in H). How does the graph H look like? 
According to our procedure every time we take x^ with the smallest index /. Therefore, 
we take the vertex %\, then connect it with itself or with the vertex 12, then, if we still have 
Xi x in (x™ 1 . . . x™* E n [g\) * we take it again and so on. At some moment we can get the 
additional (E n [g])*. This means that in this step we obtain the connected component with 
the vertexes i\, . . . , ik of the degree grater then 3. Then, using the inequality e x < l + xe x , 
we continue the procedure with (E n [g])*. At some step we can again obtain the connection 
to one of the vertexes i\, . . . ,ik, but then we must take this vertex and looking for the 
pair to it. Hence, it is easy to see that we can obtain several connected component, but 
all vertexes i p +i, . . . , i p +i lie in the same (we will call it the last) connected component of 
H thus of H . 

Since we make the finite number of steps, there are only a finite number of graphs 
H such that corresponding graphs H are equal, hence we can consider the sum over 
H instead of H. Let G {x il , . . . ,i p+ i) be the function, which corresponds to the new 
graph H . 

Note that according to (3.29) 

3+13+2 a P+l "p+l"P+2 <*P+l 

\ W* EM I* - \ W« E n [g])^ + o{l) 



n[9] 
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Therefore, we are interested in 

£ W- £a G ( Xll ,...,x tp+l ). (3.30) 

h,---,i p +i 

Since (1, . . . , 1) is an eigenvector for M* of (3.23) with eigenvalue W 2 /^, we have 

J2(M:% = —, i = -n,...,n. (3.31) 

3 

Consider any (but not the last) connected component of H . Let it has k vertexes ii, . . . , i^. 
This component is a tree with some loops and "semiloops" . The degree m^ . of each vertex 
ij is not smaller than 3, hence m := J2 m ij — 3k and the equality holds iff all degrees 
of the vertexes are 3. Let us consider the sum over ii,... ,i k - Any loop or "semiloop" 
gives W or W 1 ^ 2 respectively. Thus, since the tree has k — 1 edges, all loops give the 
contribution \Y m / 2 ~ k+1 . Using (3.31), we obtain that the contribution of the tree's edges 
is nW 2 ( fc_1 ). Therefore, since any %i j has also the coefficient W" 1 , we obtain that the sum 
over ii, . . . , i k gives n/W m / 2 ~ k+1 . Evidently m is even, and hence m/2 — k + 1 is integer. 
Moreover, since m > 3k we get m/2 — k + 1 > m/6 + 1, and thus 

m/2-£; + l>2, m < 6, 

m/2-A; + l>3, 6 < m < 12, (3.32) 
m/2-k + l>4, m > 12. 

By the same way we get that the contribution of the last connected component is not 
more than n/W m ^ 2 ~ k+1 , and m > 2{k — 1), i.e. m/2 — k + 1 > 0. In addition we have for 
the last component 

m + a>3k^k< (m + a)/2. (3.33) 

Assume that H has q connected components with the numbers of vertexes and 
with the sum of the vertex degrees bi, . . . ,b q respectively. Then (3.30) is bounded by 

1 n -r-r n 



n 



W ea W b c/ 2 ~ a i +1 J- J- vU fe j/ 2_a J +1 ' 

According to (3.32), if two of b 1: . . . , are greater than 6 or one of &i, . . . , b q ^ 1 is greater 
than 12, then the product is bounded by W~ £ ( q ~^ ■ W~ 2 , and, since b q /2 — a q + 1 > 0, 
we have that the sum (3.30) is bounded by W~ £ ( q+a \ If b q /2 — a q + 1 > 2, then we also 
have the bound W~ £<yQ+a \ Thus, it remains to consider the case b q /2 — a q + 1 = or 1 
and one of b±, . . . , b q -i is not grater then 12 and others are not grater then 6. Besides, 
bi + . . . + b q — 4ds K , since we did 2ds K steps and thus obtained 2ds K edges. Hence, since 
d > 3 and according to (3.33) b q /2 — a q + 1 < 1 yields b q < 2a, we have 

12s K < 4ds K = h + . . . + b q < 6{q - 2) + 12 + 2a = 6q + 2a. 

Therefore, q > s K or a > s K , where s K is defined in (3.28). In both cases we get the bound 
n . w~< q - 1+a ^ < n ■ W~ S£ = o(l), which gives (3.27) with (...)* and thus with (. . . )* ;0 
(see (3.29)). 

□ 
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Come back to the proof of (3.25). We can write for x G (—5, S) 

oo 

<p±(x) := exp{-ip±(x)} - 1 = 

i= 

where < (C )'. Thus 



j=—n 1=3 k=3 

where E°, S fc is the sum of all terms (Yli =1 ((f>ki x i! /W kl ))o and (n^iC^^i//^* 1 )) respec- 
tively with ki + . . . + k s — k, hi e {3, . . . , k} (E° ^ E fe) * are defined by the same way with 
(. . .)* instead of (. . .) and with \<pi\ instead of (pi). 
Denote also 

S° s = £ (fll^/^l 3 )o,*- 

ii<...<i s 1=1 

According to (4.9) (see below) we have 

s 

m^/W^) . . . (cf> ks x k : /W k *)) \ < (C ) k 5 k - 3s e Cn / w (H \x H /W\%*. 

i=i 

Hence, since the number of partitions of k to s non-zero summands is not grater than ((f) , 
we obtain 

| E 0| < eCn/W^kJ-fk^k^sgO < e Cn/W( 2Co) kJ2 5 k-3 SS 0_ (335) 

8=1 8=1 

Note now that 

— 1 2 i 4 

,o p x + px 
\x\ 3 < y - ^— , (3.36) 

and hence again according to (4.9) (see below) we get for any p > 

-1~2 /w2 i ^4 /u/4 



Besides, 



o * 

ii<...<i s i=l 



(3.37) 



1 + 9 • P ~ 1X2 2 +PX * < (1 + ^)(1 + ^) < e- 2 / 2P(1 + ^!), 

and thus, taking in account (4.9) (see below), for any p, q > such that q/p < Co with Co 
of (3.12) we have 

2n+1 ■ " - p- l x]/w 2 +pxyw^ 



fe=l j=—n 

£( ^^n (l+ ^ ))o _ 

n 4 /ti/4 



]=-n 
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where the last inequality holds in view of Lemma 3 ((. . .)o, co -q/p me ans (3.7) with 7 = 
c — q/p). This gives 

S° k < e c ^ w /q k , 
and we have from (3.37) for k > Cn/W with sufficiently big C 

S° k < e {c ™ +Cl)n/w q~ k . (3.38) 

Take q > (2|C |e) 3 . Then (3.39) and (3.38) yield for k > dn/W 

| E oj < e cn/w Y^(2C ) k 5 k - 3s q~ s < J2(S 3 q) k/3 - s < 2e Cn / w ' k . (3.39) 

s=l Q s=l 



This and (3.34) imply 



Cn/W 

(E n [<p±}) -1\<\± S*| + e~ c ^ w . (3.40) 



fc=3 



Taking into account that the number of distributions of k items into n boxes is [ n+k *) 
and using the assertion (3) of Lemma 2 we get 



I ^7 



fc / 3 ufci Jt 



/\t 1 T 3 \ 

E E X.w^^'^ 

max\xi\>SW S_1 kl <-< k ° J K---<« S 
< e -C5 2 W ^ n + ^ — 1^ < e 2felog(n/fc)-C*<5 2 U' < e -C5 2 VK/4 

where the second sum in the first line is over all collections {ki} s i=1 , = k, ki e 

{3, . . . , fc}. This yields 

S fe = S° + e~ c52w /\ ^ = ^ + e' cs2w /\ k<Cn/W 

and thus by (3.9) we have 

Cn/W Cn/W Cn/W 

1 E s °i = 1 E + e ~ C52w/A ^ E + e ' C52w/A ( 3 - 41 ) 

fc=l fe=l k=l 

Cn/W n 00 

< E K* + ^~ C52W/A < ( II (! + E - l>-,o + 2e-^ 2 ^ 4 . 

fc=l i=—n 1=3 

Since < (C ) z , there exists C such that 

00 

1 + ^10^7^ < e c ( ,3/ W 3 + ,V^) ; xe (_s\V,5W). (3.42) 
This, Lemma 3 and (3.40) yield (3.25). 
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Remark 1. 1) Using 

\fa\X < , |0 2 |=o(l), 

\ x y w < n- 1 + (n/^ 2 ) |x| 2 ^ p-V+px 4 + ^ 

instead of (3.36), we can prove (3. 41) for the series started from I = 1 with |0i| < CW~ X , 
|0 2 | = Besides, in view of (3.42) 

oo 

1 + ^ \4>l\x l /W l < e \4>i\x/W+(\4>2\-\4>i\ 2 /2)x 2 /Wi+C(x3/W3+x*/W*)^ x e (_§w,5W), 

1=1 

and hence the Cauchy-Schwarz inequality yields 

1/2 
*,0 

1/2 



( n ( x + E \<M/w l )U < ( ex p { E 2 ^(^ 3 /^ 3 + * 4 ao} 

i=—n 1=1 i=—n 

expj ]T(2|0,|x/^+ (2|0 2 | - l0i| 2 K7^ 2 })J 



i=— n 

|2 



< exp{Cm|0i| 2 + C 2 \<j) 2 \n/W}{l + o(l)) < exp{C|0 2 |n/W}, 

where to obtain the third line we use Lemma 3 for the first factor and take the Gaussian 
integral for the second factor. 

2) Define the following partial ordering. Let $i(xi, . . . , x n ), $2(^1, • • • , x n ) be two 
analytic functions in some ball with center at 0, and let the coefficients of the Taylor 
expansion of $ 2 are non-negative. Then we write 

$1 -< $ 2 (3.43) 

if the absolute value of each coefficient of the Taylor expansion of $1 does not exceed the 
corresponding coefficient o/$ 2 . 
It is easy to see that 

$3 -< $1, $4 -< $ 2 =>" $3$4 -< <&l$2- (3.44) 

The discussion above yields that if 

n 

. . . , s n ) - $ i(0, . . . , 0) -< J](l + q(si)) - 1, 

j=i 

where Sj = s(ai/W, a,i + i/W, . . . , ai + k/W, h/W, bi + i/W, . . . , bi + k/W) is a polynomial with 
s(0, . . . , 0) = 0, k is an n-independent constant, and q(s) = Ylf=i \ c jW w ^ \ c i\ — W' 1 , 
|c 2 | = o(l), |q| < (C )', / > 3, then 

\(®i)o\<(i[(l + q(smo,* + e- cs2w , 
3=1 

where s* is obtained from Sj by changing the coefficients of s to their absolute values. 
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Second step: integration over Vj 

Denote 

n 

I 



F(a,b,V) = — - V (Tr(U- n P k )*(L + A p /W)(U- n P k )i-TrLi), (3.45) 
p{ - Ao > k=-n 



dMV,A) = e J] dfi(y q ), 



Iji = J dfi(V,A). 



q=—n+l 



According to the Itsykson-Zuber formula 

JL 1 _ ( =-W 2 (a + -a_+(a 9 -6 9 )/H/)(a + -a_+K_i-6 9 _i)/PF) 

J~ = W~ An TT — = = = = . (3.46) 

9= _„+i (a+ - a- + (a, - 6g)/W)(a+ - a_ + (a,_i - 6,_i)/W) 

We want to integrate the r.h.s. of (3.21) over A). To this end we expand the 

exponent in the fourth line of (3.21) into the series of |(Vj)i 2 | 2 (note that |(V})i 2 | 2 = 
\{VM 2 i l(^)n| 2 = l(^) 22 | 2 = 1 - \{VM 2 )- K is easy to see that 

J |(y 3 ) 12 |*l) (3.47) 

1 _ e -W 2 (a+-a- + (a q ~b q )/W)(a + -a- + (a q -i^b q - 1 )/W) 



= w~ 4n yi 

q 

x(-l) 



&j (a+ - a- + (a, - &,)/W0(a+ - a_ + (a,_i - &,-i)/W0 



^ X x=W 2 (a + -a-+(a ] -i-b :j - 1 )/W)(a + -a-+(aj-bj)/W) 

If we differentiate 1 — e~ x , then we obtain e~ x which is small, since x ~ W 2 (a + — a_) 2 
(recall that |a_,-/W|, \bj/W\ < W~ K ). Therefore, we should differentiate only \jx. 

We want to prove that only the terms without |(Vj)i 2 | 2 give the contribution. Hence, 
we want to show that 



( exp{(F(a, b, V) - F(0, 0, 1))/n} - l) • III • n 2 ^ = o(l), n ->• oo, 

where III, Il 2 are the products of the Taylor's series for exp{</? + (aj/VF)} and for 
exp{ip_(bj/W)}. Besides, 



exp{(F(a,6,y)-F(0,0,y))/n}^exp{ £ ^Zw^ )' {?>M) 



C(a 3 + bj) 

3=-n 

thus we can estimate only the integral of 

exp{(F(0,0,V)-F(0,0,/))/n}. 

We can write 

e i(F(o,o,y)-F(o,o,/)) =yJ- (F(0, 0, V) - F(0, 0, J)) p . 
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Since f = J 2 + L, we have 



Tr(U. n P k )*L(U. n P k )^-TrI4 

6-6 



Tr ((U- n P k )*L (U„ n P k )L - a 2 + I 2 ) = 2a+(£ 2 - 6) • IM)i2r, 



thus 



F(0,0,V)-F(0,0,I) = 2m+ f 6) £ (|(f/_ n P fc ) 12 | 2 -|(^ t ) 12 | 2 ). (3.49) 

Pl A o) 



k=-n+l 

Hence, we have to study 



**!,..,*>,&) = (II (K^-n^-)l2| 2 - K^-n)l2| 2 ) )_, (3.50) 

where 

(...) Ji = I^J(...)dil(V,A). (3.51) 

Let p < Cn/VF for some constant C. Introduce i.i.d {tj} such that the density of the 
distribution has the form 



(a+ - a_) 2 



= g *i ex P{ - *j( a + - a -) } " io^^w^ (3.52) 



and consider the unitary matrices 



^ = ( 5' 6 1 J^i ) > (3-53) 



where 



~ = tj_ / a, - 6j \ -1/2 / a^i - 6j-i \ - V2 

J W V W(a+-a-)) \ W(a+-a_)/ 

fS- = (l-t?) 1/2 , 
and 0j G [— 7r, 7r). 

It follows from the properties of the Haar measure on the unitary group and and (3.46) 
that 



^,...,fcp(0,5):=(n(l(^-» I! ^)l2| 2 -K^-n)l2| 2 )) t 



3=1 l=-n+l 

= ® kl ,..., kp (a,b) + 0(e~ cW2 ), (3.54) 

where (. . .) t . e _ g-. means the expectation over {tj} with respect to the measure with the 
distribution (3.52) and over {9j},{9j} from — n to n. 
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Denote 



•< = 1 -( 1 + W(^K 1 + w^t))- (3 ' 55) 



Expanding Vj with respect to Sj we get 



i t 2 00 

v, = V0) + - ;)-"> - m + ^ £ v/'>«; ; 

r=l 



V?=( _ e % ), ||^||<^ (r = l,2, 



where Vj(0) is a unitary matrix (and hence ||V^(0)|| < 1), 
, . _ e ^ ^ iit7W 

and {Vj } are diagonal matrices. 

Since the integrals of e* m0J equal for m 7^ and 27r for m = 0, we conclude that if 
we replace the coefficient in front of e l6,J and e~ ldj by the bounds for their absolute values, 
then after the averaging with respect to 9j the resulting coefficients in front of s k will 
grow. Hence 



ti id. , , t 



2 



2p 



$ kl ,..,k P M ~ **!,...,*«>, 0) -< ( J] |l + JLe^Sjgfa) + ^s j9 ( Sj ) ) - 1, (3.56) 
where g(t) — Cq/{1 — Ci) with some n-independent C, Co. Moreover, 

/. 



if \ fc! 



< 



and thus we conclude 



Since for p < Cn/W we have 2p/W 2 < W~\ p 2 /W 2 = o(l), Remark 1 yields 

5 fcll ... lfcp (a,6)-$ fcli ... )fcj( (o,o)).n 1 .n 2 ' 





^2 



s (( n (1 + + w^?) - ■ n - ■ n2 -)„,. + ^ 

^ / f f ST f C P a i+}i , ^ c «T±3a \ 1 A rr tt \ , -C«5 2 w 

^ ( ( ex P { 22 [w* ■ — + ■ "wH Hl * ' Mo,* + e 

i=— n 

<-( ( -p { t {% ■ ^# + ^ • } - or ■ <t ■ <; + . 

i=— n 

where III, n 2 are the products of the Taylor's series for exp{(p + (a,j /W)} and for 
exp{(p_(bj/W)}, and n 2i * are obtained form Hi, U 2 by changing the coefficients to 
their absolute values. 
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b) - $ fcl ,..., fcp (0, 0)) • U, • n 2 ) Q < exp{|^} - 1 < exp{|^} - 1, 



We proved before that the second factor is 1 + o(l). Moreover, taking the gaussian 
integral from the first factor (similarly to Remark 1.1) we obtain that 

W 3 J ^ - L W A 
and thus, since p < Cn/W, 

Cn/W 

E ^ E n-?\(^ kl _ kp ^ 

p=l ki,...,k p 

If p >> n/VF, then 1/p! <C e -Cn/w anc i [) eilce we can change (. . .} to (. . .) 0) * (see 
Lemma 2). 

We are left to prove that 

K,.,k P (0, 0) = ( fl (l(^n)l 2 P fej (0)| 2 - |(^„)l 2 | 2 )) t = 0(1), 



where 



^(°)= n w)- 



To this end we write 

p 

H\\(U_ n P kj (0)) 12 \ 2 - \(U_ n ) 12 \ 
3=1 



l=-n+l 



< 



tj,6 



|(f/_ n P fcl (0)) 12 | 2 -|(^ n )i 2 | 



< ( |(^-„p fel -i(o)) 12 (\4 1 (o)) 22 + (^„p fcl -i(o)) 11 (Vfe 1 (o)) 12 | 2 -|(^ n ; 



12 



|(C/-nP fcl -l(0))n| 



I(^(0))l2| 



+ 



3 > v 3 >"3 
2 \(TT \ |2 



|(C/- n P fel -l(0))l2r-|(C/-„)l2 



c 



< h 

- w 2 



\(U^ n P kl ^(0)) 12 \ 2 - \(U_ n ) 12 \ 



tj ,9j ,9j 



Thus, the main term is one without |(V,-)i 2 | 2 , i.e. we can substitute Vj — 1, j — 
-n + l,...,n in (3.21). 

Hence, integrating over {V^}" = _ n+1 we obtain 

E ± = ^-8n-2 2 2n e 2(2n+l)c 0e ^« 1 -S 2 ) ^ ^ (a) /X c _ (fe) 



X 



^(2) |a j |,|fe j |<U' 1 - K 
n ra 

exp { - ¥>+(a;/W) - E M&i/W)} 



(3.57) 



j=-n 



3=-n 



i(2n+l) 



(Trt/*„Lt/_ n £-TrL£) 



(a+ - a- + (a_„ - 6_„)/W) 



x e np ( A o) 

x (a + -a- + (a n -b n )/W)dfi(U- n ) ] | rfa g d6,(l + o(l)) 



q=—n 
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Integrating over LL„ by the Itsykson-Zuber formula and using step 1 we get finally 



n 



da g d6 9 • // c+ (a) // c _(&) 



lojl.l^l^W 1 - 



q=—n 



x (a+ - a_ + (a_ n - 6_ n )/W)(a+ - a_ + (a„ - 6„)/W)(l + o(l)) 



(3.58) 



7re 2(2n+l)c (lp(Ao) 2 (47r) 2n + l g^fo _ fr)) 
^8n+2 (6 _ 6) 

By the same way we have 

^ 7re 2(2n + l) COp(Ao) 2 ( 4 7r) 2n+l sin ( n ^ 2 _ fr)) 
T ^8n+2 (6 _ 6) 

and thus (since c + = c_) 



det-'(-A + |t)|(l + o(l)) 



e 2(2n+l)co p ( Ao) 2 (47r) 2„ + l 7r 2 sin ^ _ ^ 



det 



-i 



3.3.2 E + and E_. 

Similarly to (3.21) we get 

E + = w~ 8n ~ 4 2 2n e 2( - 2n+1 ^ r++i7T ^ 1+ ^ J J 

|o J -|,|6j|<W 1 -« t/(2) 



A + 



2c+ 



(3.59) 



(3.60) 



/i c+ (a)/z c+ (6) 



x &~- 



-n+l 



(3.61) 



E (v+(a k /w)+v+(b k /w))+—i x - 1 E TV(c/_„p p )*(i p /w)(f/_„p p )| 



X 6 p——n 



x II _ t>i) 2 dn(U_ n ) Yl dii(V q )dadb(l + o(l)). 

l=—n q=—n+l 



Since 



p=—n 



exp {— E Tr (U- n P p )*(A p /W) (U- n P p )£} -< exp { ]T C(a p + & p )/nW}, 
by the same argument as for E± we get 

E = 2 2n 14/" 4 ( 2n+1 )e 2 ^ n+1 ^ r++i7r ^ 1+5 ^ 



n 



l%l,|fejl<w /1 - K 



q=—n 



x /i c+ (a)/i c+ (6) (a_ n - 6_ n ) (a n - 6 n ) (1 + o(l)) 

/n 



(3.62) 



x /i c+ (a)ii c+ (b) (a- n - &_„) (a n - 6 n ) (1 + o(l)) 
=2 2n W^" 4(2n+1) e 2{2n+1)r++i7r(5l+6) „det _1/2 £>, 
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where 



D=-A + 



2c+ 



It is easy to see (see the proof of Lemma 2) that 

\D~\ n \ = 1/deW = o(l). 

Hence, since !Rr + = Co, we get 

|E+| < 2 2n 14/~ 8ri " 4 e 2(2n+1)co+i7r(5l+ ^ ) |d et_1 - D l 
the order of E + is smaller than the order of E±. Thus we can write 



e 2( 2 n + l) C()p(Ao) 2 (47r) 2„+l 7r 2 gin ^ _ ^ 



W 8n+2 7r(6 - 6) 

This, (3.3) and (2.14) yield Theorem 1. 



(l + o(l)). (3.63) 



4 Appendix 

Proof of Lemma 1 Note that 
/*(a±) = 0, 



x=a± 



x=a± 



X 



x — 



1 - 



x 2 + A 2 /4 
1 



0, 



x=a± 



+ 



2x 2 



x 2 + A 2 /4 {x 2 + X 2 /4Y 



x=a± 



= 2(l-A 2 /4). 



Thus function f*(x) attains its minimum at a± and expanding f*(x) in x G (a± — 5, a± + 5) 
we get 

A(x) = (l-A 2 /4)( a ;-a ± ) 2 + 0(5 3 ). (4.1) 

This yields (3.13). Besides, it is easy to see, that if we take 7 = |(1 — Aq/4) then we 
obtain (3.14) for some sufficiently small S. □ 

Proof of Lemma 2 

1) Set 



-A 



1 -1 





... 









/ 


2 -1 





... 





\ 


-1 2 


-1 


... 











-1 2 


-1 


... 







-1 


2 


-1 ... 







, -A 1 = 




-1 


2 


-1 ... 







... 





-1 2 


-1 








... 





-1 2 


-1 




... 





-1 


1 


J 




V 


... 





-1 


1 


/ 



Define 



T n (x) = det (-Ai +x- 1), S n (x) = det (-A + x ■ I). 



(4.2) 
(4.3) 
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It is easy to check that 

T n (x) = (2 + x)T n _ 1 (x) -T n - 2 (x), T 1 (x) = l + x, T 2 (x) = x 2 + 3x + 1, (4.4) 
S n (x) = (1 + x)T n ^{x) - T n _ 2 (x). (4.5) 

Solving the recurrent relation (4.4) we get 

r +i + r m q , x (c m -r m )(c-i) (A R , 

T m (x) = ^—^ , S m {x) = ^—^ (4.6) 

where 

2 + x + \/x 2 + Ax 
^ ~ 2 ' 

For x = 2^/W 2 

( = 1 + ^/W + -f/W 2 + 0(W~ 3 ), W^oo. 
This and (4.5) - (4.6) yield 



T^/W 2 ) = cosh ^P(l + o(l)), S m (2l/W 2 ) = ^ sinh ^(1 + o(l)), 
and thus (3.16). Also it is easy to see that 

Gii (7) = S m (2l/W>) "727 + 

Moreover, 

{m)( W TU^M 



C~Wtanh^-^(l + o(l)) < Cl minim, W}. 
21V 



2) Take m > CW and a G R, a > 0. Note that for any sufficiently small 5 > and 

£ > 



ma,x\xi\>8W ' 

< 2^ / e J = 2 J=1 (4.7) 

i=l g=l 



e -e 2 8 2 W 2 /2 



2tt 

1=1 
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where G^ is defined in (3.17). 

Let us take e 2 = (G^ m) (a:)) _ 7 2 in (4.7). Then taking into account (3.18) and CW < 
m < 2n + 1, we obtain for a G K, a > 



7 {m) 7 {m) 



Z 



(m) 



<C 1 e~ c ^ w . 



Since m < 2n + 1, according to the first assertion of the lemma we get 



7( m ) I 

7 Ml 
J 12 I 



;i + C/Vy) m < e Cin/w , m,W ^oo. 



(4.8) 



(4.9) 



This and (4.8) yield for m > CW, 7 G C J7 > 

(m) rv(m)\ rv(rn) rv{m) rv{rri) 



zr> -z 



7 (m) 



5_£l 1 < ^3^7 . ^7 < e -C 2 S 2 W+Cn/W < e -C 3 S 2 



W 



I 



Since W 2 = n l+e we can take 5 = W~ K with k < 0/(1 + 9). 

Take now any m. Using the assertion (1) of the Lemma we can write for any e > 



(z ( m)) -i I 



-I £ ta-ay-i) 2 -^ E "' 



J=2 



J = l 



< (Z^)' 1 / e 



Xk—x\>&W 

1 / ^ i=2 ^ i=i 



x / e ^ fc + 2 

7 (fc) rv(m— k) 



I.E.fe-^-1) 2 -^ E x? m 



1 j n ^ 

g=fc+l 



< 



(m) 



3) It is easy to see that 



ax 



+ log \X\ < 1 log . 

&l 1 ~ 4 2 & a 



Thus, using assertions (1) - (2) of the Lemma and (4.9) we obtain 



|j^(m)|-l 



max |x;|>5Vy 



n 



xf ■^\x)l[dx, 

9=1 



— 1 7 



/m 



max :r;|><5W 

I7W _ 7H I 

Gifclogfc+C 2 m/H / 1 S,djy/2\ -CW& 2 

— I 7 (m) I — 

1^^7/21 

where the last inequality holds since k < Cm/W CW. □ 

Acknowledgements. I am very grateful to Thomas Spencer who drew my attention 
to this problem. 
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